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of  thermal  winds  in  the  upper  troposphere  and  lower  stratosphere  over  the  j 
White  Sands  Missile  Range  area.  A  special  series  of  about  30  concurrent 
sets  of  radiance  observations  from  the  NOAA-4  VTRR  instrument  and  wind  | 

shears  from  radiosonde  observations,  distributed  throughout  one  year,  is 
used  for  these  tests.  The  results  obtained  with  these  direct  methods  are  ; 
compared  with  results  obtained  with: 

1)  a  traditional  method,  in  which  temperature  profiles  are  first  retrieved 
from  the  satellite  radiances  and  the  thermal  winds  are  then  obtained  from 
the  horizontal  gradients  of  the  retrieved  temperatures,  and  2)  a  linear 
regression  between  observed  radiance  gradients  and  observed  wind  shears. 

The  latter  method  serves  as  an  estimate  of  the  upper  limit  of  accuracy  to 
be  obtained  by  any  method  based  on  a  linear  combination  of  radiance  gradients,  j 

The  results  indicate  that  the  direct  methods  may  be  divided  into  two  groups, 
with  much  better  retrievals  for  one  of  these  groups.  The  probable  reasons 
for  these  differences  are  identified.  The  best  direct  methods  yield  results 
comparable  to  the  traditional  method,  but  none  of  the  methods  -  not  even  the 
regression  technique  -  is  particularly  skillful.  The  lack  of  skill  in  these 
particular  cases  is  attributed  mainly  to  the  errors  associated  with  trying  j 
to  measure  relatively  small  horizontal  radiance  gradients  over  relatively 
small  horizontal  distances. 
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Foreword 

This  report  contains  a  full  description  of  our  work  on  the 
Direct  Determination  of  Wind  Shears  from  the  Gradients  of  Satellite 
Radiance  Observations  plus  reprints  of  our  paper  on  the  Impact  of 
Satellite  Temperature  Sounding  Data  on  Weather  Forecasts.  The 
latter  study  was  partially  supported  by  the  Grant. 
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ABSTRACT 


To  the  extent  that  the  wind  field  is  close  to  geostrophic, 
the  thermal  wind  is  a  good  approximation  to  the  vertical  wind  shear 
(vertical  variation  of  the  horizontal  wind).  And  since  the  thermal 
wind  is  proportional  to  the  horizontal  temperature  gradient,  the 
possibility  exists  of  determining  it  from  satellite  radiance 
observations.  Several  different  methods  are  developed  for  retrieving 
thermal  winds  directly  from  the  horizontal  gradients  of  satellite 
radiance  observations.  The  methods  are  applied  to  the  determination 
of  thermal  winds  in  the  upper  troposphere  and  lower  stratosphere  over 
the  White  Sands  Missile  Range  area.  A  special  series  of  about  30 
concurrent  sets  of  radiance  observations  from  the  NOAA-4  VTPR  instrument 
and  wind  shears  from  radiosonde  observations,  distributed  throughout  one 
year,  is  used  fo»’  these  tests.  The  results  obtained  with  these  direct 
methods  are  compared  with  results  obtained  with: 

1)  a  traditional  method,  in  which  temperature  profiles  are  first 
retrieved  from  the  satellite  radiances  and  the  thermal  winds  are  then 
obtained  from  the  horizontal  gradients  of  the  retrieved  temperatures,  and 

2)  a  linear  regression  between  observed  radiance  gradients  and  observed 
wind  shears.  The  latter  method  serves  as  an  estimate  of  the  upper  limit 
of  accuracy  to  be  obtained  by  any  method  based  on  a  linear  combination  of 
radiance  gradients. 

The  results  indicate  that  the  direct  methods  may  be  divided  into  two 
groups,  with  much  better  retrievals  for  one  of  these  groups.  The  probable 
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reasons  for  these  differences  are  identified.  The  best  direct  methods 
yield  results  comparable  to  the  traditional  method,  but  none  of  the 
.  methods  -  not  even  the  regression  technique  -  is  particularly  skillful. 

The  lack  of  skill  in  these  particular  cases  is  attributed  mainly  to  the 
errors  associated  with  trying  to  measure  relatively  small  horizontal 
radiance  gradients  over  relatively  small  horizontal  distances. 
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1.  Statement  of  the  Problem 

For  purposes  of  predicting  the  nuclear  fallout  pattern 
associated  with  the  use  of  tactical  nuclear  weapons  in  the  battlefield, 
it  is  necessary  to  know  the  vertical  distribution  of  the  horizontal 
wind  from  the  surface  of  the  earth  to  an  altitude  of  about  30  km.  In 
principle,  conventional  radiosonde  (rawinsonde)  observations  could  be 
used  for  supplying  all  of  the  required  wind  information;  however,  it  is 
not  particularly  desirable  to  use  the  radiosonde  system  for  determining 
the  wind  field  above  15  km  for  the  following  reasons: 

1)  No  other  Army  battlefield  requirement  exists  for 
meteorological  data  above  15  km.  The  nuclear 
fallout  prediction  requirement  imposes  a  substantial 
increase  (over  that  required  for  the  0-15  km 
observations)  in  meteorological  equipment,  operations 
and  manhours  on  the  battlefield. 

2)  Large  errors  are  possible  at  high  altitudes,  partly 
because  of  wind  measurement  errors  caused  by  tracking 
difficulties  at  the  low  elevation  angles  associated  with 
strong  winds,  and  partly  because  of  space  and  time 
variability  inherent  in  the  data.  Because  of  balloon 
trajectories,  the  locations  of  the  measurements  at 
altitudes  of  20  to  30  km  are  frequently  100  to  150  km 
from  the  radiosonde  balloon  launch  site. 
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To  overcome  these  problems  associated  with  the  use  of  radiosondes 
for  wind  information  above  15  km,  it  has  been  suggested^  that  a  mix  of 
conventional  radiosonde  and  satellite  radiance  observations  be  used  to 
supply  all  of  the  required  information:  the  radiosonde  would  provide  the 
wind  field  to  15  km  and  the  satellite  would  supply  the  wind  information  for 
the  15  to  30  km  layer  (which  fo^  our  purposes  we  may  call  the  stratosphere). 

The  concept  of  using  satellite  observed  radiances  to  infer  horizontal 
winds  in  the  stratosphere  is  based  on  the  following  physical  principles: 

1)  the  satellite  otserved  spectrum  of  radiance  above  a  particular 
point  on  the  earth's  surface  depends  upon  the  vertical 
temperature  profile  above  that  point. 

2)  the  vertical  variations  of  the  horizontal  wind  (windshears)  are 
directly  related  to  the  horizontal  gradients  of  atmospheric 
temperature . 

From  (1)  and  (2),  one  can,  in  principle,  estimate  the  windshears  from 
the  horizontal  gradients  of  satellite  observed  radiances. 

The  use  of  satellite  radiance  observations  for  determining  stratospheric 
wind  shears  is  particularly  attractive  due  to  the  absence  of  clouds  (which 
interfere  with  retrievals  of  temperatures)  in  the  stratosphere. 


1  Duncan,  L.D.,  1977:  SATFAL  -  The  application  of  meteorological 

satellite  data  to  nuclear  fallout  prediction.  In  Workshop  on  Satellite 
Atmospheric  Soundings  U.S.  Army  White  Sands  Proving  Ground,  SC-97. 


The  satellite  derived  product  is  a  wind  shear.  To  obtain  the  actual 
winds  in  the  stratosphere,  a  reference  level  wind  is  required:  the  15  km 
wind  observed  by  the  radiosonde  can  be  used  for  this  purpose. 

The  major  purpose  of  this  research  is  to  develop  and  test  several 
methods  for  retrieving  the  stratospheric  wind  shears  directly  from  the 
satellite  observed  radiance  gradients. 
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2.  Derivation  of  retrieval  methods. 

The  satellite  observed  radiance  for  an  observing  wavenumber 
with  negligible  contribution  from  the  earth's  surface  may  be  written 
as 


,  T  di. 


R. 

1 


i  dZ 


dZ 


(1) 


where  i  is  an  index  for  observing  wavenumber,  B  is  the  Planck  function, 
|  i  is  the  transmittance  from  the  level  Z  to  the  top  of  the  radiating 

atmosphere  Zy,  and  Z  is  a  measure  of  the  altitude  in  terms  of  scale 
heights 


Z  =  -£n(p/pQ)  (2) 


;  where  p  is  pressure  and  the  subscript  o  refers  to  the  surface  of  the 

.  earth. 

‘  > 

t-  The  horizontal  gradient  of  the  observed  radiance  in  the  meridional 

t', 

,  direction  is 

*  3 
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where  T  is  temperature  and 

dB.  di. 

1/  _ 1  _ L 

i  '  dT  dZ 


(4) 


Equation  (3)  shows  how  the  horizontal  radiance  gradient  is  related  to 
a  weighted  vertical  integral  at  the  horizontal  temperature  gradient. 


The  zonal  component  of  the  thermal  wind  in  the  layer  bounded  by 
the  levels  and  Z-j  is 

z2 

ut(Z2,Z1)  =  u(Z2)  -  u(Z1)  =  -  B.  j  (3T/3y)dZ  (5) 

^1 

where  R  is  the  gas  constant  and  f  is  the  Coriolis  plrameter. 

The  thermal  wind  is  a  good  approximation  to  the  actual  wind  difference 
between  two  atmospheric  levels  as  long  as  the  wind  field  is  dole  to 
geostrophic,  as  it  generally  is  above  the  planetary  boundary  layer. 

Comparison  of  eqs.  (3)  and  (5)  shows  that  they  both  contain 
vertical  integrals  of  the  meridional  temperature  gradient;  the 
difference  between  the  two  is  that  in  eq.  (3)  the  integral  is  over 
the  entire  atmosphere  and  is  weighted  by  K.  whereas  in  eq.  (5)  the 
integral  is  limited  to  the  layer  bounded  by  Z^  and  Z-j  and  is  unweighted. 


If  we  define  a  function  W(Z)  such  that 


W(Z)  =1  if  Z1  <  Z  <  Z2, 
W(Z)  =  0  otherwise, 


then  eq.  (5)  can  be  written 


f  utu2,Zl) 

R 


"T 


(6) 


If  eq.  (3)  is  multiplied  by  coefficients  c.  and  summed  over  the 
N  observation  wavenumbers  we  obtain 


T 

N  (  N 

I  c.r.  =  1  (E  c.K.)(3T/3y)dZ  (7) 

i=l  J0  1=1 


Inspection  of  (6)  and  (7)  indicates  that  if  we  approximate  W  by 
N 

Z  c.K.,  then  we  can  estimate  the  thermal  wind  from 
*1  1  1 

N 

ut(Z2,Z1)  =  -(R/f ) Ci r^  (8) 


The  coefficients  c^  can  be  determined  from  a  procedure  that  minimizes 
in  some  sense  the  difference  tetween  the  approximate  W-function,  which 
we  may  denote  by  W,  and  the  actual  W-function.  Fig.  1  shows  some 
actual  kernel  functions,  K.,  and  a  typical  W-function,  VI.. 


There  are  several  different  ways  in  which  the  coefficients  c^ 


can  be  determined. 


a .  Method  1 . 


One  way  is  simply  to  minimize  the  form 


rr  n  o 

|  (  E  c, K. -W j  dz 

>  i=l  1  1 
o 


This  procedure  would  minimize  the  difference  between  W  and 
W  over  the  entire  atmosphere  in  a  least  squares  sense.  The 
solution  to  this  minimization  problem  is 


C-1 

c  =  S  a 


where 


?T 

]  KiK,  dZ 

0 


K.  dZ 


(12) 


12 


With  c.j  known,  can  be  determined  from  eq.(8).  The 
meridional  component  of  the  thermal  wind,  v^,  can  be  obtained  with 
these  same  c.  by  substituting  the  zonal  radiance  gradients,  (3R^/8x), 
for  r.j  in  (8)  and  omitting  the  minus  sign. 


b.  Method  2. 


I  0 

Fleming  (1979)  (see  also  Fleming,  1972)“'  has  suggested  that 
the  coefficients  c.  be  determined  by  minimizing  the  quadratic  form 


^  \2,  N 


J  =  I  (1-W)  (  I  c . K. -W)  dZ 
i  i=l  1  1 


(13) 


subject  to  the  unit  height  constraint 


-1  [ 


2  N 


(Z  -ZJ-'  |  (  E  c  K  )  dZ 

c  1  J  i=l  1  1 

Z1 


=  1 


Fleming  refers  to  the  factor  (l-W;  as  a  penalty  function  that  tends 


1.  Fleming,  H.E.,  1979:  Determination  of  vertical  wind  shear  from  linear 
combinations  of  satellite  radiance  gradients:  A  theoretical  study.  U.S. 
Naval  Postgraduate  School  Report  NPS63-79-004,  42pp. 

2.  Fleming,  H.E.,  1972:  A  method  for  calculating  atmospheric  thicknesses 
directly  from  satellite  radiation  observations.  Preprints.  Conf.  on 
Atmospheric  Radiation,  Fort  Collins,  CO,  Amer.  Meteor.  Soc.,  134-137. 


13 


to  sharpen  up  the  approximate  W-function,  W.  Since  W=1  within 
the  layer  of  interest  and  W=0  outside  of  this  layer,  eq.  (13) 

states  that  the  integral  of  W  outside  the  layer  of  interest  should  have 
a  mininum  value.  The  unit  height  constraint  insures  that  the  integral 
of  W  is  equal  to  the  integral  of  W  within  the  layer  of  interest.  The 
solution  to  this  minimization  problem  is 

4  S  '  a 


where 


T1 

j  K.K.dZ  + 
o 


(15) 


and  i  =  Z^-Z-j  . 

c .  Method  3 


To  determine  whether  there  is  any  significant  difference  between 
the  simple  minimization  of  method  1  and  the  use  of  the  penalty  function 
as  m  method  2,  we  may  use  method  1  with  the  same  unit  height  constraint 
.is  in  method  2.  Formally,  the  solution  is  the  same  as  for  method  2, 

'  e.  eq.  (14),  but  with  S.^  defined  as  in  method  1,  i.e.  eq.  (11). 
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d.  Method  4. 


Methods  2  and  3  botn  use  the  unit  height  constraint.  This 
insures  that  areas  such  as  those  labelled  C  in  figure  2  would  cancel 
each  other.  However,  the  excess  areas  labelled  A  and  B  are 

2 

undesirable  and  lead  to  inaccuracies  in  the  solution.  Fleming  (1972) 
suggests  compensating  for  these  areas  by  multiplying  the  solution  of 
method  2  by  the  factor 


f 


1/  WdZ 

j 


o 


06) 


e.  Method  5. 


To  insure  that  the  correct  solution  is  obtained  in  the  case 
where  the  thermal  wind  is  not  a  function  of  altitude,  i.e.,the 
horizontal  temperature  gradient  is  constant  with  altitude,  we  may 
proceed  as  follows.  Combine  (5),  (6)  and  (8)  to  obtain 


h  n 

|  OT/ay)dZ  *  E 


c.r. 
1  i 


]t  N 

i  (  1 

>  i  =  l 
o 


ciK|(9T/3y)dZ 


07) 


where  the  caret,  as  before,  indicates  the  estimated  value.  For 
(3T/9y)  constant  with  altitude. 
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OT/3y)(Z2-Z1)  =  (3T/.1y) 


(18) 


From  (IS)  it  can  be  seen  that  the  estimated  temperature  gradient 
3T/3y,  will  equal  the  actual  temperature  gradient  if 


(z2-zi) 


-1 


(19) 


This  is  a  different  form  of  the  unit  height  constraint.  Eq.  (19)  can  be 
interpreted  as  stating  that  the  vertical  integral  of  the  approximate 
W-function  should  be  equal  to  the  vertical  integral  of  the  exact 
W-function.  This  seems  to  be  a  more  appropriate  normalization  constraint 
than  the  one  given  in  method  2.  The  use  of  this  constraint  with  the 
minimization  of  eq.  (9)  constitutes  method  5.  The  solution  is 

,  l-  bT  S"1  a  , 

c  *  S' '  a  +  — — v— —  5  b  (20) 

b  S'  b 


where 


bi 


K.dZ 


(21) 


-!  %l  4  a.  . 
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f.  Traditional  and  regression  methods. 

The  above  five  methods  may  be  called  direct  retrieval 
methods  since  the  thermal  winds  are  obtained  directly  from  the  radiance 
gradients.  Using  a  set  of  colocated  radiance  and  radiosonde  data,  we 
shall  compare  the  results  obtained  from  the  different  direct  methods. 

We  shall  also  compare  the  results  of  the  direct  methods  with  results 
obtained  from  method  T,  the  traditional  method,  in  which  the  thermal 
wind  is  computed  from  the  temperatures  retrieved  from  radiance  obser¬ 
vations,  and  from  a  linear  regression  method  based  on  the  relationship 
between  the  observed  wind  shears  and  the  observed  radiance  gradients. 
The  results  obtained  from  the  regression  method  represent  a  limiting 
accuracy  for  a  particular  data  set  for  any  direct  linear  method. 
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3.  Experimental  details. 

A  special  series  of  31  concurrent  (within  1  hour)  sets  of 
radiance  observations  from  the  NOAA-4  VTPR  instrument  and  wind 
shear  from  radiosonde  observations,  distributed  throughout  one 
year  (Feb.  to  Dec.  75),  for  the  White  Sands  Missile  Range  area 
in  New  Mexico  (see  Fig.  3  is  used  for  evaluating  the  different 
methods.  The  horizontal  radiance  gradients  over  this  area  are 
estimated  from  a  7  x  7  slightly  rectangular  array  of  satellite 
radiance  observations  extending  over  an  area  of  about  500  km  on 
the  side.  The  spot  size  (50%  power  contour)  of  each  radiance 
observation  ranges  from  60  km  at  the  sub-satellite  point  to  about 
80  km  at  the  end  of  each  scan  line.  To  compute  the  radiance 
gradients,  four  averages  of  3x7  grid  points  each,  as  shown  in 
Fig.  3,  were  made  of  the  radiances  and  of  the  location  coordinates. 
The  radiance  gradients  are  then  obtained  from  application  of  the 
Taylor  approximation  formulas  to  these  average  values 

ARi  =  W  If  Axl  +  If  ^  (22) 

aR2  =  Rc  '  Rd  =  If  Ax2  +  ly  Ay2 


(23) 


where  R^,  Rg,  Rg  and  R_.  are  the  average  radiances  (at  a  particular 
wavenumber)  for  the  regions  A,  B,  C  and  D,  as  shown  in  Figs.  3  and 
4,  and  Ax-j ,  Ay-j ,  AX2  and  are  shown  geometrically  in  Fig.  4. 

These  two  simultaneous  equations  can  be  readily  solved  for  the 
unknown  radiance  gradients  3R/8x  and  9R/9y.  This  simple  method 
gave  approximately  the  same  results  as  a  planar  fit  by  a  regression 
technique. 

"Ground  truth"  was  obtained  from  radiosonde  balloons  launched 
(nearly)  simultaneously  from  the  three  stations  shown  in  Fig.  3 
about  one  hour  prior  to  the  overpass  of  the  NOAA-4  satellite.  The 
average  of  the  three  measured  wind  profiles  was  used  as  "truth"  data 
for  the  different  methods. 

It  was  decided  to  compare  the  direct  methods  using  just  the 
three  upper  atmosphere  channels  (695.5,  679.5,  667.8  cm”  ) .  Sample 
comparisons  with  the  use  of  four  channels  indicated  no  improvement  in 
the  results. 

The  atmospheric  transmission  for  each  of  CO^  band  observing  channel 
is  assumed  to  be  a  product  of  the  CO^ »  ^0  and  03  transmissions  in  that 
channel.  A  zenith  angle  correction  to  the  transmissions  is  also 
included.  These  transmissions  are  computed  for  each  case  from  a 
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first  guess  temperature  and  water  vapor  profile  and  from  a 
standard  ozone  distribution. 

In  the  traditional  method,  temperature  retrievals  are 
obtained  with  an  inversion  scheme  similar  to  that  of  Hogan  and 
Grossman  (1972). 


I 


I 

31 
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4.  Results 

Before  discussing  the  results,  it  is  instructive  to  compare 
the  approximate  W-functions  of  the  different  direct  methods.  Fig.  5 

A 

shows  the  W-function  for  the  five  direct  methods  for  the  layer 

A 

125-25  mb.  The  W-functions  of  methods  2  and  3  are  almost  identical. 
This  suggests  that  there  is  no  particular  advantage  to  be  gained  by 
using  the  penalty  function  as  in  eq.  (13)  rather  than  a  straightforward 
minimization  technique.  It  will  be  recalled  that  methods  2  and  3 
both  use  the  unit  height  constraint  applied  to  the  layer  of  interest. 

A 

While  this  insures  that,  within  the  layer  of  interest,  the  area  under  W 
is  equal  to  the  area  under  W,  it  also  causes  greater  side  lobes  outside 
the  layer  of  interest  than  those  of  the  other  methods.  Although 
method  4  also  makes  use  of  the  unit  height  constraint  in  the  layer  of 
interest,  the  final  solution  coefficients  are  reduced  by  a  common 
factor  in  order  to  reduce  the  side  lobes,  thus  nullifying  the  intent  of 

A 

of  the  unit  height  constraint.  The  W  functions  of  the  other  layers 
have  similar  characteristics. 

Just  as  the  W-functions  of  the  different  methods  group  themselves 
into  two  distinct  groups  -  those  with  unit  height  constraint  within 
the  layer  of  interest  and  those  without  (included  here  is  method  4)  - 


so  do  the  results.  There  were  no  substantial  differences  in  the 
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results  obtained  with  methods  2  and  3,  and  there  were  no  substantial 
differences  in  the  results  obtained  with  methods  1,  4  and  5.  There¬ 
fore,  in  the  discussion  that  follows  the  results  of  method  1  are  used 
to  represent  methods  1,  4  and  5,  which  we  may  call  group  1,  and  the 
results  of  method  2  are  used  to  represent  methods  2  and  3,  which  we 
may  call  group  2. 

Table  1  compares  the  arithmetic  (or  systematic)  errors  of  the 
group  1  and  group  2  direct  methods  and  the  traditional  (T)  method. 
Also  shown  in  the  table  are  the  observed  mean  wind  shears.  The 
wind  shears  represent  the  change  in  wind  speed  from  the  base  to  the 
top  of  each  layer. 

In  general,  the  group  2  methods  have  a  much  greater  systematic 
error  than  the  group  1  methods,  especially  for  the  u-component  of 
the  wind  shear.  The  arithmetic  errors  of  group  1  are,  in  general, 
somewhat  less  than  those  of  the  traditional  method.  The  arithmetic 
errors  of  group  1  and  the  traditional  method  are  small  compared  to 
the  observed  mean  u-component  wind  shears,  but  are  of  the  same  size 
as  the  observed  v-component  shears,  which  are  of  the  order  of  1  m  s  ^ 
for  each  layer. 

Table  2  shows  a  comparison  of  the  root-mean-square  errors  of 
the  various  methods.  Also  included  in  this  table  are  the  standard 
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errors  (SI.)  of  a  linear  'egression  method  based  on  the  relationship 
between  the  observed  wind  shears  and  observed  radiance  gradients. 

These  standard  errors  represent  estimates  of  the  minimum  root-mean-square- 
errors  that  can  be  achieved  by  any  direct  linear  retrieval  method.  In 
addition,  the  table  includes  the  standard  deviations  (a)  of  the  observed 
wind  shears 


inspection  of  Table  2  reveals  that  the  root-mean-square-errors 
of  group  2  are  typically  more  than  twice  as  large  as  those  of  group  1, 
the  mean  root-mean-square-errors  for  the  four  layers  between  250  and 
10  mb  oei  nq  2i  m  s  '  and  3ms  resnecti  ve  ly  *  for  the  u-compcfient ,  and 
7  5  and  15  5ms  '  for  the  v-component 

A  comparison  of  group  1  root-mean-square-errors  with  those  of 
the  trad’t'ional  method  shows  that  they  are  about  the  same  magnitude  for 
u-component  out  greater  for  group  1  for  the  v-component.  However,  for 
neither  method  does  tne  root-mean-square- error  go  much  below  the 
standard  deviation  of  the  observations.  This  implies  that  the  error 
obtained  by  assuming  that  the  wind  shear  is  always  equal  to  the  observed 
c1 imatolog’cal  value  at  White  Sands  would  be  no  greater  than  that 
obtained  when  satellite  observations  are  used  to  estimate  the  wind 


shea*  s . 


Comparison  of  the  root-mean-square-errors  of  the  direct 
methods  and  the  standard  errors  of  the  regression  technique  shows 
that  the  difference  between  the  two  is  larger  for  the  v-component 
of  the  wind  shear  than  fo>'  the  u-component . 

In  both  Tables  1  and  2,  a  separate  row  of  results  is  presented 
for  the  125-25  mb  layer.  This  layer  is  approximately  twice  as 
large  as  the  other  layers,  and  it  was  thought  that,  in  view  of  the 
"resolving  power"  difficulties  of  the  basic  radiance  observations, 
the  results  might  be  improved  for  such  a  "thick"  layer.  However, 
inspection  of  the  tables  indicates  that  there  is  no  particular 
reduction  in  the  errors  for  this  thicker  layer. 
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5.  Conclusions  and  discussion. 

Several  methods  have  been  derived  for  determination  of 
thermal  winds  directly  from  the  horizontal  gradients  of  satellite 
radiance  observations.  Tnese  methods  are  intercompared  using 
a  special  series  of  satellite  VTPR  radiance  and  radiosonde  wind 
observations  over  the  White  Sands  Proving  Ground  area  in  the  United 
States.  The  direct  methods  are  also  compared  to  the  "traditional" 
method,  in  which  thermal  winds  are  computed  from  the  horizontal 
temperature  gradients  obtained  from  temperature  profiles  retrieved 
from  the  satellite  radiance  observations,  and  to  a  linear  regression 
method 

Of  the  direct  methods  tested,  those  that  do  not  apply  the  unit 
height  constraint  to  the  layer  of  interest  give  the  best  results.  The 

A 

use  of  the  unit  height  constraint  causes  larger  side  lobes  in  the  W- 
function,  thus  introducing  more  undesirable  information  from  the 
region  outside  the  layer  of  interest  than  in  the  case  of  the  methods 
that  have  no  such  constraint. 

There  appears  to  be  no  particular  advantage  to  introducing  a 
penalty  function  (1-W)  ,  into  the  integral  form  that  is  minimized 
to  obtain  the  solution  coefficients  for  the  direct  methods. 


The  best  direct  methods  gave  results  comparable  to  those 
obtained  with  the  traditional  method.  However,  all  the  methods 
examined  are  not  particularly  skillful  at  retrieving  the  wind 
shears  from  the  radiance  observations.  This  lack  of  skill  is 
probably  associated  with  a  number  of  factors. 

We  believe  that  the  major  factor  leading  to  lack  of  skill  is 
the  error  associated  with  trying  to  evaluate  a  relatively  small 
horizontal  gradient  over  a  relatively  small  horizontal  distance. 
For  the  geographical  region  and  atmospheric  layers  studied  in 
this  research,  typical  u-component  wind  shears  are  of  the  order 
of  10  m  s'1 .  Such  wind  shears  imply,  through  the  thermal  wind 
equation,  meridional  temperature  gradients  of  about  3.5°C  per 
1000  km.  Differentiation  of  the  Planck  function  with  respect  to 
temperature  indicates  that,  for  the  temperatures  and  wavenumbers 
of  concern  here,  a  change  of  1°C  is  equivalent  to  a  change  of  one 
radiance  unit  (mW/m  sr  cm  ).  Thus,  the  observed  wind  shears 
are  associated  with  radiance  gradients  that, are  of  the  order  of 
3.5  radiance  units  for  1000  km.  However,  in  the  present  investi 
gation  these  radiance  gradients  are  estimated  from  radiance 
observations  over  characteristic  horizontal  distances  of  about 
300  km  -  which  means  radiance  differences  of  about  1  radiance 
unit  or  temperature  differences  of  about  l°r  must  be  measured. 


26 


The  \l TPR  instruments  had  characteri Stic  sensitivities  of  0.25 
radiance  units  (even  larger  at  667.8  cm  ^),  a  significant 
fraction  of  the  radiance  differences  that  had  to  be  measured  in 
order  to  retrieve  the  observed  thermal  winds.  Thus,  the  VTPR 
instruments  are  just  not  sensitive  enough  to  measure  the  small 
horizontal  radiance  differences  over  distances  of  ~  300  km 
associated  with  typical  meridional  temperature  gradients  in  the 
stratosphere  over  White  Sands.  The  situation  Tor  the  v-component 
is  even  worse,  since  typical  v-component  wind  shears  are  less  than 
those  of  the  u-component. 

Other  factors  contributing  to  errors  in  the  retrievals  include 
the  differences  between  the  actual  and  the  approximate  W-functions, 
uncerta inties  in  transmission  functions,  and  possible  differences 
between  the  thermal  wind  shear  and  the  actual  wind  shear. 

The  above  discussion  leads  us  to  suggest  that  much  better 
results  might  be  obtained  by  applying  the  direct  methods  of  the 
group  1  type  to  regions  of  the  atmosphere  with  typically  greater 
wind  shears  than  those  encountered  in  the  stratosphere  over  White 
Sands  and/or  to  characteristic  horizontal  distances  greater  than 
the  300  km  used  in  the  present  study.  Improved  instrumentation 
in  the  current  generation  of  satellite  sounding  equipment  should  also 
lead  to  improved  retrievals. 
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Figure  Captions 


Fiq.  1  K.  and  W  functions. 

°  i 

Fig  2  Schematic  W  and  K  functions. 

Fig.  3  Typical  array  of  satellite  radiance  observation  points 
over  the  White  Sands  area.  Each  dot  represents  the 
center  of  a  spot  observed  by  the  satellite.  The  circled 
crosses  show  the  locations  of  the  radiosonde  stations  that 
were  used  for  "ground  truth".  The  rectangles  A,B,C  and  D 
are  used  in  the  computation  of  radiance  gradients,  as 
explained  in  the  text  (see  also  fig.  4). 

F 1 g .  4  Schematic  diagram  showing  how  observations  from  areas  A,G,C 
and  D  are  used  to  compute  radiance  gradients.  See  text. 

/N 

Fig  5  Comparison  of  approximate  (W)  and  actual  W  functions  for 


125-25  mb  layer. 
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Table  1  Mean  observed  wind  shears  and  mean  arithmetic  errors  of 
wind  shear  retrievals  ("'  s'^). 

Method  1  represents  tht  group  of  methods  1,  4  and  5; 
Method  2  represents  the  group  of  methods  2  and  3;  and 
N  is  the  number  of  compari sons . 


Layer 

u- 

component 

j 

V 

-component 

1 

(mt) 

(km) 

Observed 

Method 

1  2 

1 

T 

Method 

Observed  1  2 

- ^ 

i 

T 

N 

250-100 

i 0  4-16.2  : 

-8.2 

6.9  7.7 

2.4* 

i 

-0.4 

0.7 

1  .0 

0.0* 

31 

100-50 

16.2-20  6 

1 

-11.6 

1.9  -7  .4 

5.8 

-1.0 

0.5 

0.4 

1.1 

31 

50-20 

20.6-26.5  1 

1  -4.0 

-1  .7  -9.8 

-1.6 

-0.3 

-0.2 

-1  .0 

0  7 

28 

20-10 

26.5-31 . 1 

1.2 

-0.2  4.1 

-3.4 

1.1 

-0.7 

-0.9 

-0.7 

15 

125-25 

14.8-25.0 

-19.8 

3.4  -1.9 

7.9 

-2.1 

1  .3 

1  0 

2.3 

31 

*  Values  for  Method  T  are  for  the  300-100  mb  layer  (9  6-16.2  km). 
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Table  2.  Root-mean-square-errors ,  standard  error  of  regression  technique  (S.E.), 
and  standard  deviation  of  observations  of  wind  shear  (c)  (m  s  ^ ) . 

Method  1  represents  the  group  of  methods  1,  4  and  5; 

Method  2  represents  the  group  of  methods  2  and  3,  and  N  is  the  number 
of  comparisons. 


L.aye” 


u-component 


250-100  10.4-16.2 


14  7.6* 


v-component 


11.  14  8.3* 


100-50  16.2-20.6 


14  7.9 


6.5  11  3.5 


30-20  20.6-26.5 


22  5.2 


8.0  18  4.4 


20-10  26.5-31.1 


35  6.3 


4.5  19  2.5 


8.9  21  6.7  6.5  7.9  8.4  19  4.7  4.4  6.0 


125-25  14.8-25.0  j  8.7  10  10.1  6.4  9.4  I  12.  15  5.5  5.5  8.0  1  31 

*  Values  for  Method  T  are  for  the  300-100  mb  layer  (9.6-16.2  km) 


MS 
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5  Comparison  of  approxima*,*-  (V)  anl  ai-tnal  V  f-irotlors  for 
l?5-?5  mb  layer. 
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Abstract 

The  concept  of  radiometric  sounding  of  atmospheric  tempera¬ 
ture  profiles  from  satellites  was  first  demonstrated  with  data 
gathered  by  infrared  spectrometers  on  the  Nimbus-3  satellite 
in  1969.  Operational  satellite  sounding  over  oceanic  areas  was 
introduced  by  the  VTPR  (Vertical  Temperature  Profile 
Radiometer)  instrument  on  the  NOAA  2  satellite  in  1972. 
Karly  evaluations  of  these  new  observational  data  centered  on 
their  accuracy  compared  to  data  obtained  from  the  con¬ 
ventional  radiosonde  system.  More  recent  evaluations  have 
focused  on  the  impact  of  the  satellite  temperature  soundings 
on  numerical  weather  forecasts.  In  this  paper,  we  review  the 
results  of  such  impact  tests  in  several  countries.  On  the 
average,  the  inclusion  of  satellite  sounding  data  leads  to  a 
small  improvement  in  the  numerical  forecasts. 


1.  Introduction 

The  VTPR  instrument  on  the  operational  l  .S.  weather 
satellites  provided  about  1000  sounding  temperature 
profiles  per  day  over  the  oceans  of  the  globe.  This 
number  is  about  one  half  of  the  total  number  of  radio¬ 
sonde  reports  available  daily.  The  VTPR  temperature 
profiles,  processed  in  W  ashington,  D.C.,  were  available 
over  the  global  meteorological  telecommunications 
circuit.  The  infrared  and  microwave  sounding  systems 
on  the  L.S.  experimental  satellite  Xiinl>tis-6  (prototypes 
of  the  sounding  system  for  the  operational  satellite 
TIROS-N),  were  capable  of  providing  about  5000  tem¬ 
perature  soundings  daily,  over  land  and  water.  Satellite 
temperature  soundings  are  in  the  form  of  temperature 
versus  pressure  and  can  he  converted  to  thicknesses 
directly  or  to  geopotential  heights  with  the  use  of  the 
hypsometric  equation  and  a  knowledge  of  the  height  of 
some  reference  level,  generally  the  1000  mb  surface. 
Satellite  data  are  asynoptic ;  the)  can  lie  included  in  the 
analysis  for  a  particular  observing  time  either  by  means 
of  an  observing  window  (e.g.,  the  use  of  data  ±3  hours 


1  Expanded  version  of  paper  presented  at  the  COSPAR 
Symposium  on  the  Remote  Sounding  of  the  Atmosphere  from 
Space,  Innsbruck,  Austria,  1-3  June,  1978. 
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around  the  observing  time)  or  an  asynoptic  assimilation 
scheme.  Most  evaluations  of  the  impact  of  satellite 
soundings  on  forecasts  are  made  by  running  two  parallel 
sets  of  forecasts,  one  set  based  on  initial  analyses  that 
do  not  include  satellite  temperature  sounding  data,  the 
other  on  initial  analyses  that  include  these  data.  The 
two  sets  of  forecasts  are  then  verified  against  the  ob¬ 
served  states  of  the  atmosphere,  generally  in  regions 
with  dense  radiosonde  networks. 

A  number  of  such  impact  tests  have  been  performed 
by  various  groups  around  the  world.  The  results  have 
been  published  in  scientific  journals  or  in  reports.  We 
have  attempted  to  collect  the  results  of  these  studies  for 
presentation  here.  As  part  of  this  attempt,  we  have 
written  to  a  number  of  national  meteorological  services 
requesting  any  information  they  could  provide  on  the 
impact  of  satellite  soundings  on  forecasts.  Part  of  the 
information  provided  was  qualitative  in  nature.  We 
summarize  in  the  next  section  both  the  quantitative  and 
qualitative  information  we  have  compiled. 

2.  Results 

a.  Quantitative  impart 

Impact  evaluation  tests  have  been  performed  by  a 
number  of  groups,  using  different  analysis,  forecasting, 
and  verification  methods.  We  have  chosen  the  48  h 
500  mb  geopotential  height  forecast  for  evaluating  the 
impact  of  satellite  sounding  data.  This  particular  fore¬ 
cast  was  chosen  because  it  represents  an  important 
atmospheric  level  and  forecast  range,  and  results  for 
this  forecast  were  available  from  almost  all  of  the  im¬ 
pact  studies.  The  most  popular  forecast  verification 
criteria  were  the  rms  errors  and  SI  skill  score.  We  shall 
present  the  impact  results  in  terms  of  one  or  tile  other 
(or  both,  where  available)  of  these  criteria.  The  -SI 
score  is  given  by 

100  L  Ion  |  EIC.il 

where  eu  is  the  error  of  the  forecast  height  difference 
and  (>l  is  the  observed  or  forecast  height  difference  be¬ 
tween  grid  points,  whichever  is  larger.  The  lower  the 
score,  the  more  skilful  the  forecast.  Twenty  years  of 
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Tablk  1.  Impact  of  satellite  soundings  on  rms  errors  of  48  h  500  mb  geopotential  height  forecasts  (m). 
(Positive  impact  represents  reduction  in  rms  height  eiror.) 


Number  of 


Source 

Verification  area 

Season 

Data 

forecasts 

NOSAT 

SAT 

Impact 

Desmarais  et  al.  (1978) 

Eastern  N.  America 

Summer 

V  +  N 

10 

45.8 

47.9 

-2.1 

Western  N.  America 

43.5 

4S.0 

-1.5 

X.  America 

Winter 

V  +  N 

15 

65.0 

63.7 

1.3 

Halem  el  al.  (1978) 

N.  America 

Winter 

V  -1-  N 

11 

77.9 

72.8 

5.1 

Bonner  et  al.  (1976) 

N.  America 

Spring 

V 

9 

63.4 

60.6 

2.8 

Atkins  and  Junes  (1975) 

Europe,  Atlantic, 

Spring 

V 

7 

94.0 

88.0 

6.0 

Canada,  XE  USA 

Druyan  et  al.  (1978) 

Europe,  Mid- East, 

X.  Africa 

Experimental  VTPR 

Winter 

V 

26 

70.0 

67.0 

3.0 

Operational  VTPR 

Winter 

V 

13 

69.3 

67.3 

2.0 

Kelly  (1977)* 

Australia 

Winter 

V 

9 

49.0 

43.6 

5.4 

V  =  VTPK  soundings,  N  =  Nimbus-6  soundings. 
*  36  h  forecasts. 


experience  in  the  use  of  SI  scores  at  the  l  .S.  National 
Meteorological  Center  (NMC)  indicates  that  Un¬ 
practical  range  between  a  perfect  and  worthless  forecast 
at  500  nth  is  20  70  (Fawcett,  1077). 

Table  1  shows  the  impact  of  satellite  sounding  data 
on  numerical  forecasts  as  measured  by  the  reduction  in 
rms  geopotential  height  errors,  a  positive  impact  repre¬ 
senting  a  reduction  in  the  error  when  satellite  sounding 
data  are  included.  NOSAT  refers  to  the  rms  error  when 
no  satellite  soundings  are  included  in  the  initial  analyses 
and  SAT  refers  to  those  cases  in  which  satellite  sound¬ 
ings  were  added  to  the  data  base  of  the  initial  analyses. 
Kxcept  for  the  slimmer  results  of  Desmarais  et  al.  (1078) 
the  impacts  are  all  positive,  indicating  improved  fore¬ 
casts  when  satellite  sounding  data  are  included.  How¬ 
ever,  the  impacts  are  all  small,  ranging  from  3  6  in. 
while  the  error  levels  of  the  forecasts  are  in  the  range  of 
45  100  m.  The  largest  improvement,  percentagewise, 
is  obtained  in  the  Southern  Hemisphere  (Kelly,  1077) 
where  the  impact  is  greater  than  10r,’ .  In  all  studies,  the 
small  average  impact  obtained  is  the  result  of  a  combina¬ 
tion  of  positively  influenced  and  negatively  influenced 
forecasts  and  not  the  result  of  a  consistent  impact  on 
each  forecast. 


Table  2  shows  impact  results  in  terms  of  SI  skill 
scores,  a  positive  impact  meaning  a  reduction  in  the  SI 
score  when  satellite  sounding  data  are  included:  a  re¬ 
duction  in  SI  represents  an  improvement  in  the  fore¬ 
cast.  The  SI  results  confirm  the  rms  error  results.  Small 
positive  impacts  are  obtained  in  all  studies,  except  for 
the  tiny  negative  impact  in  the  summer  results  of 
I  lesmarais  et  al.  (1978).  The  most  impress/  ve  results  are 
once  again  in  a  Southern  Hemisphere  study  where  a 
reduction  of  almost  5  points  is  registered  in  the  SI 
score  when  satellite  soundings  are  included  (Kelly  et  al.. 
1978).  Once  again  each  impact  listed  consists  of  the 
mean  of  a  set  of  forecasts,  some  of  which  showed  positive 
impact  and  some  of  which  showed  negative  impact. 

It  is  of  interest  to  compare  the  winter  results  of 
Desmarais  et  al.  (1978)  at  NMC  with  those  of  Halem 
et  al.  (1978)  at  HISS  (Goddard  Institute  for  Space 
Studies).  Ten  of  the  forecasts  were  common  to  both 
studies.  The  NMC  obtained  a  positive  impact  of  0.5 
points  in  SI  score,  representing  a  reduction  of  SI  from 
34.8  to  34.3.  C.ISS  obtained  a  positive  impact  of  1.9 
points  in  SI.  representing  a  reduction  from  39.6  to  37.7. 
Halem  et  al.  (1978)  imply  that  C.ISS’s  greater  positive 
i  in  pact  is  the  result  of  a  better  assimilation  and  analysis 


Tablk  2.  Impact  of  satellite  soundings  on  SI  scores  of  48  h  geo|x>tential  height  forecasts. 
(Positive  impact  represents  reduction  in  SI,  improvement  in  forecast.) 


Source 

Verification 

area 

Season 

Data 

Number  of 
forecasts 

NOSAT 

SAT 

Impact 

Desmarais  et  al.  (1978) 

X.  America 

Summer 

V  +  N 

10 

46.0 

46.1 

-0.1 

Winter 

V  +  \ 

15 

34.8 

34.3 

0.5 

Halem  et  al.  (1978) 

X.  America 

Winter 

V  +  N 

11 

39.6 

37.7 

1.9 

Bonner  et  al.  (1976) 

N.  America 

Spring 

V 

9 

44.2 

45.6 

1.4 

Kelly  (1977)* 

Australia 

Winter 

V 

9 

42.0 

40.0 

2.0 

Kelly  et  al.  (1978)** 

Australia 

Winter 

N 

28 

42.2 

37.4 

4.8 

V  =  VTPR  soundings,  N  =  Nimbus  6  soundings. 
*  36  h  forecasts. 

**  24  h  forecasts. 
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f'  to.  1.  Mean  annual  SI  skill  score  of  NMC’s  36  h  forecasts 
at  500  mb  during  the  period  1955  76.  (From  Halem  el  at., 
1978:  extension  of  data  presentcrl  by  Fawcett .  1977.) 


Fortcast  p«riod  (h ) 


Kig.  2.  Si  skill  score  of  NMC  six-level  primitive  equation 
model  at  500  mb  as  a  function  of  forecast  period.  (After 
Desmarais  et  at.,  1978.) 


scheme.  Desmarais  et  at.  (1978)  suggest  that  the 
greater  impact  registered  by  < » ISS  is  simply  due  to  more 
room  for  improvement  in  C.ISS’s  analysis  and  forecast 
system,  as  evidenced  by  better  (by  5  points)  NOSAT  SI 
scores  of  XMCs  forecasts.  (That  the  difference  in  SI 
scores  is  not  due  to  different  areas  or  grid  intervals  in  the 
two  studies  is  shown  by  Desmarais  et  at.  (1978)  who  used 
exactly  the  same  verification  code  for  evaluating  (iISS 


and  NMC  SI  scores,  and  found  the  same  ~5  SI  point 
difference.)  I  n  particular.  Desmarais  et  al.  (1978)  imply 
that  CISS’s  analyses  without  satellite  soundings  are 
poorer  than  NMC’s  and,  therefore,  the  addition  of 
satellite  soundings  can  improve  (.ISS’s  analyses  more 
than  they  can  improve  NMC’s  analyses;  hence,  a 
greater  impact  in  C.ISS’s  forecasts.  To  help  resolve  this 
question,  it  would  be  desirable  to  apply  C.ISS’s  assimila¬ 
tion  and  analysis  scheme  to  the  N  MC  forecast  model  in  a 
series  of  impact  tests,  experiments  along  these  lines  are 
currently  underway  (t.hil  et  at.,  1979a). 

The  question  naturally  arises  as  to  the  meaning  or 
significance  of  a  small  improvement,  say  1  or  2  points, 
in  the  SI  skill  score  at  500  mb.  Some  light  can  be  shed 
on  this  question  by  Figs.  1  and  2.  Figure  1  shows  the 
mean  annual  SI  score  of  the  NMC’s  56  h  forecasts  at 
500  mb  during  the  period  1955  76  (from  Halem  et  til., 
1978 ;  extension  of  data  presented  by  Fawcett,  1977).  A 
general  trend  toward  lower  SI  scores,  i.e.,  improved 
forecasts,  can  be  seen  during  the  period.  Over  the  last 
10  years  there  has  been  about  a  lOCj  improvement  in 
the  SI  score,  i.e.,  a  reduction  of  about  4  points.  If  this 
rate  can  also  be  assumed  representative  of  the  48  h 
forecast,  then  a  reduction  of  1  or  2  points  in  the  SI 
score  represents  an  improvement  that  currently  takes 
about  2  to  5  years  to  achieve.  Figure  2  (after  Desmarais 
et  til..  1978)  shows  the  500  mb  SI  score  of  the  NMC 
6-level  primitive  equation  model  as  a  function  of 
forecast  period.  For  forecast  periods  less  than  60  h.  the 
forecast  deteriorates  at  a  rate  of  about  5  SI  points  for 
12  h.  Thus,  an  improvement  of  1  or  2  points  in  the  SI 
score  is  equivalent  to  an  extension  of  forecast  capability 
by  about  2}  5  li  at  the  48  h  forecast  range. 

A  1  or  2  point  SI  impact  may  be  small  compared  to 
the  differences  in  the  performance  of  differing  analysis 
and  forecasting  systems,  e.g.,  (>ISS  versus  NMC. 
Operational  experience  at  NMC  (Shuman,  1978)  has 
shown  step  jumps  greater  than  1  or  2  points  when  im¬ 
proved  models  are  adopted  (e.g.,  alter  introduction  of 
the  6-level  primitive  equation  model,  or  the  fine  mesh 
model).  But  improvements  in  the  forecast  models  may 
be  more  difficult  to  achieve  today  than  in  the  past  (there 
is  less  room  for  improvement)  and  a  1  or  2  point  im¬ 
provement  in  SI  score  may  be  significant  at  today’s 
state  of  model  development. 

In  all  the  studies,  as  previously  noted,  the  small 
average  impact  obtained  is  the  result  of  a  combination 
of  positively  influenced  and  negatively  influenced  fore¬ 
casts  and  not  the  result  of  a  consistent  impact  on  each 
forecast.  This  is  illustrated  (Kig.  3)  by  a  breakdown  of 
the  impacts  registered  by  Halem  eta/.  (1978)  for  North 
America  during  a  winter  period.  Although  the  average 
impact  was  positive,  representing  a  reduction  of  1.9 
points  in  skill  score,  in  3  out  of  the  11  forecasts  the 
impact  was  negative. 

b.  Qualitative  impact 

Qualitative  information  on  the  impact  of  satellite 
soundings  on  numerical  weather  forecasts  was  obtained 
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Impact  on  SI 

Fig.  3.  Histogram  of  impacts  of  satellite  sounding  data  on 
48  It  500  mb  SI  scores  for  North  America  during  a  winter 
period.  (After  Halem  etal.,  1078.) 


from  several  national  meteorological  services.  Their 
findings  are  summarized  below. 

Weather  Bureau.  Republic  of  .South  Africa  The  1000  30i 
mb  thickness  Held  is  one  of  two  primary  fields  predicted, 
and  both  satellite  sounding  data  and  satellite-derived 
winds  are  included  in  the  operational  manual  analyses 
of  this  field.  The  Weather  Bureau.  Republic  of  South 
Africa,  indicates  that  "although  we  have,  as  yet.  made 
no  objective  assessment  of  the  impact  of  satellite  sound¬ 
ing  data  and  cloud  imagery  on  the  ipiahty  of  our 
numerical  analyses  and  numerical  prognoses,  we  are 
convinced  that  it  must  lie  considerable”  ( Triegaardt. 
1978.  private  communication). 

Meteorological  .Service.  France.  The  experience  of  the 
French  Meteorological  Service  indicates  that  if  one 
removes  the  satellite  sounding  data  from  the  initial 
conditions  of  a  particular  forecast,  the  resulting  forecast 
is  not  significantly  different  from  a  forecast  based  on 
initial  conditions  including  the  satellite  data  (Mittner. 
1978,  private  communication).  However,  if  one  re¬ 
moves  the  satellite  soundings  from  the  operational 
objective  analysis  scheme  for  a  period  of  several  con¬ 
secutive  days,  one  observes  a  progressive  degradation 
of  the  analyses  and  the  forecasts.  particular!'  in  those 
areas  not  well  covered  by  other  types  of  observations. 

Meteorological  .Service,  Federal  Republic  of  (iermany.  The 
FRC.  Meteorological  Service  indicates  that  satellite 
sounding  data  are  used  routinely  in  the  numerical 
analysis  procedure.  Although  no  impart  experiments 


have  been  performed,  "we  feel  that  the  old  system  with 
eight  weather  ships  in  the  Atlantic  Ocean  provided 
more  valuable  information  about  actual  weather  phe¬ 
nomena  which  cannot  be  replaced  by  currently  avail¬ 
able  satellite  soundings"  (Buschner,  1978,  private 
commit  nication). 

3.  Discussion  and  conclusions 

Our  survey  of  the  impact  of  satellite  soundings  on 
in  merical  forecasts  indicates  that  current  satellite  data 
as  used  in  current  analysis  and  forecast  models  produce, 
on  average,  a  small  improvement  in  the  numerical  fore¬ 
casts.  (It  should  be  noted  that  the  \imbus-6soundings 
were  never  used  operationally  lint  only  in  simulations 
of  operational  forecasts).  However,  this  small  average 
positive  impact  is  not  based  on  consistent  small  positive 
impacts  in  each  forecast,  but  rather  on  an  average  of 
forecasts  with  positive  impact,  negative  impact  and  no 
impact.  The  small  effect  of  the  satellite  data  is  also 
evidenced  by  almost  all  qualitative  and  synoptic  evalua¬ 
tions  of  their  impact.  Subjective  evaluations  of  NOSAT 
and  SAT  numerical  forecasts  by  experienced  fore¬ 
casters  (Desmarais  et  a!..  1978)  indicated  little,  if  any. 
difference  between  the  two  for  the  summer  test  period; 
during  the  winter  test  period,  in  those  cases  wlcre  the 
SAT  and  NOSAT  forecasts  differed  (only  26%  of  ail 
the  cases),  the  SAT  forecast  was  judged  the  better  fore¬ 
cast  aimed  j  of  the  time.  From  subjective  evaluation  of 
the  usefulness  and  accuracy  of  prognostic  charts  gen¬ 
erated  from  SAT  and  NOSAT  forecasts,  and  from 
verifications  of  precipitation  forecasts  as  might  be 
issued  in  local  forecast  operations  based  on  outputs  front 
SAT  and  NOSAT  forecasts.  Halem  et  al.  (1978)  find 
"beneficial  but  modest  impacts."  However,  Halem  et  al. 
(1978)  also  find  that  in  those  regions  of  the  forecast  area 
where  the  impacts  are  relative!'  large  (>  |96  mi  in 
500  mb  height  for  72  h  forecasts),  the  impact  is  positive 
in  seven  out  of  nine  cases.  In  Atkins  and  Jones  (1975), 
in  which  positive  impact  was  registered  by  the  satellite 
soundings,  “on  no  occasion  would  the  forecast  of 
weather  issued  for  the  British  Isles  have  been  affected 
by  the  inclusion  of  satellite  lata.”  Druyatt  et  al.  (1978). 
in  discussing  a  forecast  for  which  the  mis  height  errors 
were  significantly  reduced  by  inclusion  of  VTI’R  ob¬ 
servations  in  the  initial  analysis,  state  that  “we  would 
hi  si  talc  to  claim  that  the  satellite  data  changed  a  poor 
forecast  into  an  accurate  one."  The  pit  ture  is  not  as  bad 
in  the  Southern  Hemisphere.  Here,  the  few  available 
studies  suggest  that  satellite  soundings  produce  greater 
improvements  in  the  numerical  forecasts. 

Possible  reasons  for  lack  of  impact  have  been  dis¬ 
cussed  by  Desmarais  et  al.  (1978).  Tracton  and 
McPherson  (1977).  and  Drux an  etal.  (1978).  Satellite 
soundings  show  differences  of  1  5°t'  when  compared  to 
radiosonde  observations.  Temperature  analyses  that 
include  satellite  data  show  smaller  variance  than 
analyses  based  on  radiosonde  data  alone.  Figure  4  (after 
Desmarais  et  al.  (1978))  shows  tin  consistently  smaller 
;s)d>  available  potential  merge  (a  measure  of  tempera- 
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Etc.  4.  Eddy  available-  potential  energy  —  AK  (101  Jm  *) 
of  the  atmosphere  between  8S0  and  200  mb,  from  20-90°N  for 
SAT  and  NOSAT  analyses.  (After  Desmarais  et  at.,  1978.) 

ture  variance)  of  analyses  that  include  satellite  sound¬ 
ings.  The  satellite  sounding  data  must  compete  with  an 
initial  analysis  that  is  derived  from  other  sources.  Ap¬ 
parently  the  quality  of  Northern  Hemisphere  initial 
analyses,  even  over  data-sparse  regions  such  as  ocean s, 
is  not  improved  substantially  by  the  addition  of  the 
satellite  soundings  (inridently.  what  is  usually  thought 
of  as  a  data-sparse  region  may  have  significant  inputs 
from  aeroplane  observations,  winds  derived  from  satel¬ 
lite  observations,  and  a  reasonably  good  12  It  forecast) ; 
hence,  there  is  a  lack  of  impact  on  the  forecasts.  In  the 
Southern  Hemisphere,  where  the  quality  of  the  initial 
analyses  is  not  nearly  as  good  as  that  in  the  Northern 
Hemisphere,  the  satellite  data,  despite  problems  of 
quality,  lead  to  improved  analyses  and  forecasts. 

The  above  discussion  takes  a  rather  pessimistic  view 
of  the  impact  of  satellite  soundings,  at  least  in  the 
Northern  Hemisphere.  The  optimists  may  argue,  how¬ 
ever,  that  the  small,  average  improvement  obtained  is 
significant  when  compared  to  the  rate  at  which  fore¬ 
casts  have  improved  over  the  past  decade,  that  it  is 
becoming  more  and  more  difficult  to  improve  the  ac¬ 
curacy  of  forecasts,  and  that  even  a  small  improvement 
may  be  significant.  They  can  argue  that  satellite  sound¬ 
ings  have  led  to  improvements  in  the  forecasts  that 
might  have  taken  2  5  years  to  achieve  otherwise,  and 
that  forecast  capability  has  been  extended  by  2}  5  h 
at  the  48  h  range.  They  may  have  a  point. 

One  should  not  take  the  present  conclusions  as  the 
final  word  on  the  subject.  The  number  of  impact  studies 
thus  far  conducted  is  relatively  small,  and  the  number 
of  rases  in  each  study  has  been  barely  enough  to  provide 
a  general  picture  of  impacts.  More  impart  studies  are 
required.  Satellite  soundings  are  not  point  observations 
like  radiosonde  observations,  but  volumetric  observa¬ 
tions:  satellite  soundings  have  their  own  error  charac¬ 
teristics.  The  weather  services  of  the  world  have  not 
had  much  experience  with  the  use  of  this  new  type  of 
data.  Improved  analysis  schemes  (see.  for  example. 
Phillips.  1976)  may-  enhance  the  impart  of  satellite 
soundings.  There  is  also  the  possibility  of  reducing  the 


errors  of  satellite  soundings  (through  increased  knowl¬ 
edge  of  atmospheric  transmission  characteristics,  for 
example),  although  the  basic  problem  of  poor'  vertical 
resolution  will  always  remain  for  passive  observing 
techniques.  Smith  (1977,  private  communication)  in¬ 
dicates  that  the  satellite  observations,  which  have  a 
horizontal  field  of  view  of  the  order  of  50  km,  have  been 
processed  to  yield  a  temperature  sounding  representa¬ 
tive  of  a  horizontal  area  of  the  order  of  a  500  km  square, 
Asa  result  of  this  "horizontal  averaging,"  the  horizontal 
fields  of  the  satellite  data  are  smoother  and  contain  less 
structure  than  those  typically  obtained  from  analysis 
of  radiosonde  observations.  Smith  (1977,  private  com¬ 
munication)  believes  that  processing  the  satellite  data 
at  their  —  50  km  resolution  with  a  feyy  hundred  kilo¬ 
meters  separation  between  processed  soundings  would 
yield  horizontal  fields  more  nearly  compatible  with 
radiosonde  observations  and  with  the  numerical  weather 
prediction  models  as  currently  formulated.  However,  by 
virtue  of  the  inherent  vertical  averaging  of  the  satellite 
temperature-sounding  technique,  a  certain  amount  of 
horizontal  smoothing  is  automatically  introduced  into 
the  analyses.  In  addition,  the  statistical  climatological 
nature  of  the  retrieval  techniques  may  also  suppress 
horizontal  structure.  To  what  extent  the  processing  of 
the  satellite  data  at  their  basic  ~  50  km  horizontal 
resolution  can  restore  the  actual  horizontal  structure  of 
the  temperature  field  remains  to  be  seen.  In  any  case, 
the  TOVS  (TIROS  Operational  Vertical  Sounder)  on 
the  third  generation  polar  satellite  system  of  NOAA. 
the  TIROS  N  (the  first  TIROS-N  was  launched  on  18 
October  1978).  should  provide  a  wealth  of  data  for 
further  impact  studies. 

Note  added  in  proof:  At  the  time  this  pat>er  was  yvritlcn,  the 
results  of  the  recent  American  impact  tests  were  available  only- 
in  report  form.  Summaries  of  these  studies  are  beginning  to 
ap,._  r  in  tbi*  scientific  journals.  Some  of  the  GISS  results 
appear  in  t  <)/.  (I979h1.  One  new  result  is  based  upon 

subjective  cr a. nation  by  experienced  forecasters  of  NOSAT 
and  SAT  forecasts  for  the  winter  test  period.  The  GISS  fore¬ 
casters  find  that,  in  those  cases  yvhere  there  were  differences 
in  the  SAT  and  NOSAT  forecasts  of  500  mb  geopotential 
height  (46','.  of  total  cases!,  the  SAT  forecast  was  the  better 
one  about  J  of  the  time.  The  NMC  impact  work  is  reported 
Tracton  et  nl.  (1979a)  and  Tract  on  et  ol.  (1979b). 

TIROS-N,  launched  on  15  October  1978,  yvith  the  third 
generation  of  vertical  temperature  sounding  instruments  on 
board,  is  capable  of  providing  10  000  soundings  a  day.  How¬ 
ever,  since  redundant  observations  in  polar  regions  are  not 
processed,  the  total  is  reduced  to  8000  |ier  day,  of  which  2000 
are  transmitted  via  the  Global  Telecommunications  System 
(GTS).  The  TIROS-N  vertical  sounding  system  is  described 
by  Smith  et  nl.  ( 1979),  and  an  early  evaluation  of  the  TIROS-N 
temperature  soundings  is  given  by  I’hillips  et  nl,  (1979).  Both 
papers  appear  in  this  issue  of  the  lUl.I.l-TIN .  Impact  studies 
are  planned  for  the  future. 

TIROS-N  was  joined  in  s|iare  on  27  June  1979  by  N’OAA-6, 
the  second  in  this  new  series  of  l\  S.  0|>crutional  orbiters. 
When  the  NOAA-6  sounding  system  becomes  operational,  the 
two  satellites  will  provide  4000  soundings  per  day  over  the 
GTS. 

The  GISS  group  that  performed  the  NASA  impact  studies 
is  now  at  the  Laboratory  for  Atmospheric  Sciences,  NASA 
Goddard  Space  Flight  (  enter. 
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